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Abstract 

We present the calculations of the complete next-to-leading order (NLO) QCD effects on the 
single top productions induced by model-independent tqg flavor-changing neutral-current couplings 
at hadron colliders. Our results show that, for the teg coupling the NLO QCD corrections can 
enhance the total cross sections by about 60% and 30%, and for the tug coupling by about 50% and 
20% at the Tevatron and LHC, respectively, which means that the NLO corrections can increase 
the experimental sensitivity to the FCNC couplings by about 10%— 30%. Moreover, the NLO 
corrections reduce the dependence of the total cross sections on the renormalization or factorization 
scale significantly, which lead to increased confidence on the theoretical predictions. Besides, we 
also evaluate the NLO corrections to several important kinematic distributions, and find that 
for most of them the NLO corrections are almost the same and do not change the shape of the 
distributions. 

PACS numbers: 14.65.Ha, 12.38.Bx, 12.60.Cn 



'Electronic address: csli@pku.edu.cn 



I. INTRODUCTION 



The top quark is the heaviest particle so far discovered, with a mass close to the elec- 
troweak (EW) symmetry breaking scale. Thus it is a wonderful probe for the EW breaking 
mechanism and new physics beyond the standard model (SM) through its decays and pro- 
ductions at colliders. One of the most important aspects of the top quark physics is the 
investigation of the anomalous flavor-changing neutral-currents (FCNC) couplings in the top 
quark sector. Within the SM, the FCNC couplings are absent at the tree level, and occur 
at high order in perturbation theory through loop diagrams which are further suppressed 
by the Glashow-Iliopoulos-Maiani (GIM) mechanism 1|. On the other hand, these FCNC 
couplings can be enhanced to observable levels in some new physics models [2|, such as the 
models with extra quarks 3J, two Higgs doublet models [4|, supersymmetric models [5], ex- 
tra dimensions models 6], little Higgs models Q, or technicolor models As the coming 
Large Hadron Collider (LHC) will produce abundant top quark events (about 10 8 per year), 
even in the initial low luminosity run (~ 10 fb _1 /year) 8 x 10 6 top quark pairs and 3 x 10 6 
single top quarks will be produced yearly, one may anticipate the discovery of the first hint 
of new physics by observing the FCNC couplings in the top quark sector. 

Since we do not know which type of new physics will be responsible for a future deviation 
from the SM predictions, it is necessary to study the top quark FCNC processes in a model- 
independent way by an effective Lagrangian. In general, any new physics at a high scale 
A will manifest themselves at energies below A through small deviations from the SM, 
which can be described by an effective Lagrangian containing higher dimensional SM gauge 
invariant operators jjj. As we know, the dimension five operators break baryon and lepton 
number conservation, the lowest order operators considered are dimension six. For the tqg 
anomalous FCNC couplings, the only independent one is the chromo-magnetic operator 
induced by a dimension six gauge invariant operator before spontaneous gauge symmetry 
breaking [101 ] . 

9s E ^i^uT a qG% + H.c, (1) 



q=u,c 



where T a are the Gell-Mann matrices, G a are the field strength tensors of the gluon, and 
K t q (<? = u i c ) are rea l coefficients that define the strength of the couplings. To test the above 



couplings at the Tevatron and 



single top quark production 
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iC, the direct top quark production 



ll| and the FCNC 



13] are the most promising processes, and the top quark 



rare decay process is not so efficient compared with the above two processes due to the 
large SM backgrounds [14J. Currently the most stringent experimental constraints for the 
tqg anomalous couplings are /cf u /A < 0.018 TeV -1 and ftf c /A < 0.069 TeV -1 , given by the 
CDF collaboration through the measurements of the direct top quark production [jjj] using 
2.2 fb _1 of data. The DO collaboration also analyzes 230 pb -1 of data and provides similar 
constraints, K 9 tu /K < 0.037 TeV -1 and /sf c /A < 0.15 TeV -1 , based on the measurements of 



the FCNC single top production 



16| using the next-to-leading order (NLO) K factors of the 



direct top quark production process 17] . Taking into account the difference of the amount 



of data, we can see that the discovery potential of these two channels are almost the same at 
the Tevatron. As for the LHC, according to the tree-level analysis in Ref. 13], for the FCNC 
single top quark production, kL/A can be detected down to 0.0061 TeV -1 and 0.013 TeV -1 
with an integrated luminosity of 10 fb -1 for q = u and c, respectively, thus it is one of the 
most important channels to detect the top FCNC couplings at the LHC. 

As we know, the leading order (LO) cross sections for processes at hadron colliders suffer 
from large uncertainties due to the arbitrary choice of the renormalization scale (jj, r ) and 
factorization scale (p.f), thus are not sufficient for the extraction of the FCNC couplings. 
Then high order corrections including NLO QCD and resummation effects are needed in 
order to improve the theoretical predictions. The NLO QCD and threshold resummation 
effects for the direct top quark production are studied in Ref. [l^], and the NLO QCD 
corrections to the top quark rare decays via FCNC couplings can be found in Ref. 18]. 
As for the FCNC single top production, all the predictions are still at the LO. In this 
paper, we present the complete NLO QCD corrections to the single top quark production 
via model-independent tqg FCNC couplings at both the Tevatron and LHC. 

The arrangement of this paper is as follows. In Sec. [TT]we show the LO results. In Sec. IHIt 
we present the details of the NLO calculations, including the virtual and real corrections. 
Section IIVI contains the numerical results, and Section |V] is a brief summary. 

II. LEADING ORDER RESULTS 

At the LO there are three main subprocesses which contribute to the single top production 
via the FCNC couplings at hardron colliders: 

g q >tg, g g >tq, q(q, q') q ► t q(q, q'), (2) 



where q is either u quark or c quark. The corresponding Feynman diagrams are shown in 
Fig. [TJ It should be noted that, after expanding the effective operator in Eq. ([I]), there are 
two kinds of FCNC vertices, one is the three point vertex, another is the four point vertex, 
and both of them contribute to the subprocesses at the LO. 
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FIG. 1: The LO Feynman diagrams for the single top quark production via the FCNC couplings. 
The ones for gg and qq subprocesses can be obtained by crossing symmetry. 



The LO squared amplitudes for the corresponding subprocesses in four dimensions are 

A 2 



\M*\ 2 gq (s,t } u) 
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\M^ gq (s,t,u) 



\M^\ 2 qlq (s } t } u) 



+ (8s 3 + 92ts 2 + 92t 2 s + 8t 3 ) m 10 - 2 (4s 4 + 53ts 3 + 109t 2 s 
+53t 3 s + 4t 4 ) m 8 + (4s 5 + 66ts 4 + 234t 2 s 3 + 234t 3 s 2 + 66t 4 s 
+4t 5 ) m 6 - st (19s 4 + 122ts 3 + 221t 2 s 2 + 122t 3 s + 19t 4 ) m 4 

+6s 2 t 2 (5s 3 + 14ts 2 + 14t 2 s + 5t 3 ) m 2 - 3s V (5s 2 + Ats + 5t 2 )) 

-4A 2 

(3m 6 - 3(3s + 2t)m 4 + 2 (6s 2 + 8ts + 3t 2 ) m 2 - 2s (3s 2 

+5ts + 5t 2 )) , 
-4A 2 
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(m 2 (s + t) - 2st) 



(3) 



"g \M*\ 2 gg (u,t,s), \M*\ 2 qq (s,t,u) = \MZ\ 2 qq (t,s,u), 

respectively, where m is the top quark mass, and A = 87rcv s /«L/A, the colors and spins of the 
outgoing particles have been summed over, and the colors and spins of the incoming ones 
have been averaged over, s, t, and u are Mandelstam variables, which are defined as 

8 = (pi + p 2 ) 2 , t = (px- p 3 ) 2 , u = (pi - p 4 ) 2 - (4) 



After the phase space integration, the LO partonic cross sections are given by 

The LO total cross section at hadron colliders is obtained by convoluting the partonic cross 
section with the parton distribution functions (PDFs) Gj/p for the proton (antiproton): 

o B = ^2 / dx ^ dx 2 [ G a/p 1 {xi, fif)G b/ p 2 {x 2 , fi f )a^ b ] , (6) 

ab J 

where ///is the factorization scale. 



III. NEXT-TO-LEADING ORDER QCD CORRECTIONS 

The NLO corrections to the single top production via the FCNC couplings consist of the 
virtual corrections, generated by loop diagrams of colored particles, and the real corrections 
with the radiation of a real gluon or a massless (anti) quark. We carried out all the cal- 
culations in the 't Hooft-Feynman gauge and used the dimensional regularization (DREG) 



scheme 191 in n = 4 — 2e dimensions to regularize all the divergences. Moreover, for the real 

□ 

corrections, we used the dipole subtraction method with massive partons [20] to separate 
the infrared (IR) divergences, which is convenient for the case of massive Feynman diagrams 
and provides better numerical accuracy. 



A. Virtual corrections 



The virtual corrections for the single top production via the FCNC couplings include the 
box diagrams, triangle diagrams and self-energy diagrams as shown in Figs. [2]l3] and Fig. H] 
for the gq and qq initial state subprocess, respectively. For simplicity we did not show the 
diagrams that only differ by the exchange of identical external particles. The loop diagrams 
for the gg and qq initial state subprocess can be obtained from Figs. [2}|3] and Fig. H] by 
crossing symmetry, and the ones for q'q are just part of Fig. HI 

All the ultraviolet (UV) divergences appearing in the loop diagrams are renormal- 
ized by introducing counterterms for the wave functions and mass of the external fields 
(5z!f\ 8Z^\ 8Z2 , 8m) , and the coupling constants (8Z gs ,8Z K a / A ). We define these coun- 
terterms according to the following conventions. For the external fields, we fix all the 
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FIG. 2: One-loop Feynman diagrams for the subprocess gq — > if?, part I. 



renormalization constants using on-shell subtraction, and, therefore, they also have IR sin- 
gularities: 



2tt c V 3 2J \e uv e IR J 6tt e e uv ' 
a, / 1 1 



3vr \ e c/v e /R 
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37T \e[/y 

where C e = T(l + e)[(47r/i^)/m 2 ] e and n/ = 5. For the renormalization of g s , we use the MS 
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FIG. 3: One-loop Feynman diagrams for the subprocess gq — > tg, part II. The big gray circle 
represents the gluon self-energy diagrams. 



scheme modified to decouple the top quark 



211 ]. i.e. the first rif light flavors are subtracted 



using the MS scheme, while the divergences associated with the top quark loop are subtracted 
at zero momentum: 

Z J e uv LZn e uv 



5Z n 



^r(i + £)( 4.).U 



(8) 



Thus, the renormalized strong coupling constant a s evolves with rif light flavors in this 
scheme. Finally, for the renormalization constants of the FCNC couplings SZ K 9 / A , we adopt 



the MS scheme and adjust it to cancel the remaining UV divergences exactly: 



a. 



5Z </A =£r(l + 6)(47T) 
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FIG. 4: One-loop Feynman diagrams for the subprocess qq — » 



and the running of the FCNC couplings are given by 22] 

_ <(m) /a s (m)V /(3/3o) 
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a a (ji r ) 

with /3b = 11 - 2n//3. 

The squared amplitudes of the virtual corrections are 



M 



' gq(q q )\l-loop 



^22Re(M l °°P> i M B *) gg{gq) + 2Re(M con M B *] 



(10) 



(11) 



where M l °°^ denote the amplitudes for the i-th loop diagram in Figs. [2JI3] or Fig. HJ arid 



M 



gq r qq \ are the corresponding counterterms. All the one-loop integrals in t 



le loop amplitudes 



can be calculated using the standard Passarino-Veltman techniques [23|, and the explicit 
expressions for the scalar integrals containing IR divergences can be found in Ref. 24|. In 
Eq. (ITT)) , all the UV divergences cancel each other, leaving the remaining IR divergences 
and the finite terms. Because of the limited space, we do not shown the lengthy explicit 
expressions of the virtual corrections here. The IR divergence of the virtual corrections can 



S 



be written as 

\M\ 2 ag \i-ioo P ,iR = -^-D e \ x \M B \ 2 + terms proportional to — 1 , 
3tt [e IR e IR ) 

\Wji-ioo P ,iR = -—DA-%- x \M B \\ Q + terms proportional to — 1, (12) 

7T {e IR e IR ) 

where D e = [(47r/i^)/s] e /r(l — e), and |M S | 2 are the squared Born amplitudes given in 
Eq. Q. The 1 /em terms can not be factorized in a trivial way due to the nontrivial color 
structures of the LO amplitudes, and can be expressed as combinations of the LO color 
correlated squared amplitudes as shown in the next section. 

B. Real corrections 

At the NLO the real corrections consist of the radiations of an additional gluon or massless 
(anti) quark in the final states, including the subprocesses 

g q ► t g g, t q q, t q' <f , g g ► t q g, q(q, q') q >t q(q, q') g, (13) 

as shown in Fig. [5j It should be noted that in our NLO calculations we did not include the 
contributions from the SM on-shell production of the top pair with subsequent rare decay 
of one top quark, pp(p) — > ti — > t + q + g, which provide the same signature as the single 
top production via the FCNC couplings and can be calculated separately. We will discuss 
these contributions in the numerical results. 

Before performing the numerical calculations, we need to extract the IR divergences in 
the real corrections. In the dipole formalism this is done by subtracting some dipole terms 
from the real corrections to cancel the singularities and large logarithms exactly, and then 
the real corrections become integrable in four dimensions. On the other hand, these dipole 
subtraction terms are analytically integrable in n dimensions over one-parton subspaces, 
which give e poles that represent the soft and collinear divergences. Then we can add them 
to the virtual corrections to cancel the e poles, and ensure the virtual corrections are also 
integrable in four dimensions. This whole procedure can be illustrated by the formula [20I ]: 



& NLO 



I [{d*«)^-(da A ) e J + I \da v + Ida 

J m+l Jm L J\ 



(14) 



e=0 



where m is the number of final state particles at the LO, and da A is a sum of the dipole 
terms. Besides, at hadron colliders, we have to include the well-known collinear subtraction 
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FIG. 5: Typical Feynman diagrams of the real corrections, others not shown can be obtained by 
the exchange of the identical external particles. 



counterterms in order to cancel the collinear divergences arising from the splitting processes 
of the initial state massless partons. Here we use the MS scheme and the corresponding 
NLO PDFs. 

For the process with two initial state hadrons, the dipole terms can be classified into four 
groups, the final-state emitter and final-state spectator type, 



Aj,fc(Pl, — ,Pm+l) 
1 



(Pi+Pj) 2 -mtj 
the final-state emitter and initial-state spectator type, 



' i J 



m,a\-- • i 



ij, ...;a, 



(pi +Pj) 2 - rn\ Xij )a 
the initial-state emitter and final-state spectator type, 

Vf (pi, ...,p m+1 ;p a , ...) = 

1 1_ _ ~ ~. ,T,T 



2 j I ■••> J ; ■ ) u,t > ■■■Im.aii 



(15) 



X T ~ _ 

T 2 3 V?-|...,ij,...;o,...)m,a, (16) 



(17) 
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and the initial-state emitter and initial-state spectator type, 



V ai '\ Pl ,...,p m+1 - Pa ,p b ) 



-,(...; ai,b\ 



T h T 



ai V^\...;ai,b) -j 



(18) 



where a, 6 and i, j, ... are the initial and final state partons, and T and V are the color charge 
operators and dipole functions acting on the LO amplitudes, respectively. The explicit 
expressions for x i)abl Xy j0 and V can be found in Ref. |20J. The integrated dipole functions 
together with the collinear counterterms can be written in the following factorized form 



m,ab 



J d§ {m \p a ,p b ) m ,ab{-;Pa,Pb\1m+a+b{^)\-;Pa,Pb)i 
+ y dx d& m) (x Pa , Pb ) 

m,a' &(...; xp a , pb\P'^ L +b {x) + K^_ 6 (x)|...; Xp a) Pb)m,a' 

a' J ° J 

+{a «-> b), 



(19) 



where x is the momentum fraction of the splitting parton, d^ m ^ contains all the factors 
apart from the squared amplitudes, I, P, and K are insertion operators defined in [20]. 
For simplicity, in all the above formulas we do not show the jet functions that define the 
observables and are included in our numerical calculations. 

The operators P and K provide finite contributions to the NLO corrections, and only 
the operator I contains the IR divergences 
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(20) 



with 



i fQ 2 ^ e 



V(sj k ,mj,mk;e IR ) = — 

"Ujk \ Sjk 

7j 



1-& 1-2^ 1 



IR 



c 

Fj(0,6 IR ) = — , ^(mj 7^ 0,e JK ) = — , 



(21) 



where Cf = 4/3, 7 9 = 2, and 7 9 = 11/2 — n//3. And Sjfc, fjfc, and p n are kinematic 



11 



variables defined as follows 



s^-fc = 2pjp fc , Q, fc = s jk + m- + m fc , v jk = \ 1 - 



2 ; 



l-^ + 2m^-mj _^ 

l + ^ + 2m^-mj-mD ln ^ 

When inserting Eq. (120]) into the LO amplitudes for the gq and gg subprocesses as shown 
in Eq. (Tl9l) . we can see that the IR divergences, including the 1/epR terms, can be written 
as combinations of the LO color correlated squared amplitudes and all the IR divergences 
from the virtual corrections in Eq. (TT2l) are canceled exactly, as we expected. 

IV. NUMERICAL RESULTS 

In the numerical calculations, we investigate the NLO QCD effects on the total cross 
sections, the scale dependence, and several important distributions at both the Tevatron 
and LHC For the single top production via the FCNC couplings at the NLO, the final state 
consists of a top quark plus one or two partons which may form jets. We use the cone jet 
algorithm with AR = 0.5, and require the jet to have pt > 20GeV and \r)\ < 2.5. We 
require that there is at least one jet in the final state, which means the signal we considered 



is t+jet + X, unlike the direct top production. All the input parameters are taken to be 25]: 



m t 



171.2TeV, a s {M z ) = 0.118, k 9 J A = k 9 JA = O.OlTeV" 1 . (23) 



The running QCD coupling constant is evaluated at the three-loop order 



3 



and the 



CTEQ6M PDF set [26( is used throughout the calculations of the NLO (LO) cross sec- 
tions. Both the renormalization and factorization scales are fixed to the top quark mass 
unless specified. We have performed two independent calculations for the virtual correc- 
tions and the integrated dipole terms, and used the modified MadDipole [27j package to 
generate the Fortran code for the real corrections. The numerical results of the two groups 
are in good agreement within the expected accuracy of our numerical program. 
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FCNC coupling tug (LO) tug (NLO) teg (LO) teg (NLO) 



LHC 



re/A 



■O.OlTeV- 



pb 



6.77 



,41 



1.10 



1.49 



Tevatron ( q^jT^t ) 2 fb 86 



129 



6.2 



10.2 



TABLE I: The LO and NLO total cross sections for the single top quark production via the FCNC 
couplings at both the LHC and Tevatron. 
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FIG. 6: The LO total cross sections and NLO K factors as functions of the leading jet transverse 
momentum cut. 



In Table HI we list some typical numerical results of the LO and NLO total cross sections 
for the single top quark production via the FCNC couplings. We can see that the cross 
sections can reach about 8 pb and 1 pb at the LHC, for the tug and teg coupling, respectively. 
But if we take the FCNC coupling values to the experimental upper limits 15j, then the 
cross sections can reach as large as several tens of pb. In Fig. [6] we present the LO total cross 
sections and the K factors unlo/^lo as functions of the leading jet transverse momentum 
cut, respectively. It can be seen that, for the teg coupling the NLO corrections can enhance 
the total cross sections by about 60% and 30%, and for the tug coupling by about 50% and 
20% at the Tevatron and LHC, respectively. And the K factors decrease with the increasing 
transverse momentum cut. In Sec. IIIIBI we mentioned that the SM on shell production of the 
top pair with subsequent decay also contributes to the same final state as the FCNC single 
top production. According to Ref. fl8|. the former contributions are about 7 x 10~ 4 pb and 
8 x 10~ 2 pb at the Tevatron and LHC, respectively, for our chosen FCNC coupling values. 
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Compared with the results in Table HI we can see that these contributions are negligibly 
small for the tug coupling, but can reach about 10% of the LO total cross sections for the 
teg coupling at both the Tevatron and LHC. 
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FIG. 7: Scale dependence of the total cross sections at the Tevatron, the black lines represent the 
LO results, while the red ones represent the NLO results. 




l-i/m, ^m t 

FIG. 8: Scale dependence of the total cross sections at the LHC, the black lines represent the LO 
results, while the red ones represent the NLO results. 

In Figs. [7]and [8]we show the scale dependence of the LO and NLO total cross section for 
three cases: (1) the renormalization scale dependence fi r = = m t , (2) the factorization 

scale dependence /i r = m t , fj,f = /jl, and (3) total scale dependence \i r = = fi. It can 
be seen that the NLO corrections reduce the scale dependence significantly for all three 
cases, which make the theoretical predictions more reliable. For example, at the Tevatron 
for the tug coupling, when the scale \i varies from 0.2m t to 5m t , the variations of the total 
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cross sections are 90% and 40% for case (1), 50% and 20% for case (2), and 150% and 
50% for case (3) at the LO and NLO, respectively. Note that, from Fig. [8] it seems that 
the NLO corrections did not reduce the factorization scale dependence at the LHC This 
is because the PDFs of the two incoming partons have opposite trends with the increasing 
factorization scale and the scale dependence cancel each other at the LO, but there is no 
such cancellation at the NLO. In order to further explain this scenario, in Fig. [9] we plot 
the factorization scale dependence of the gu initial state subprocess at the LHC, which 
contributes about 90% of the LO total cross section for the tug coupling. It can be seen 
that, if we fix the factorization scale of one incoming parton and only change another one, 
the NLO corrections indeed improve the factorization scale dependence. 
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FIG. 9: Factorization scale dependence for the gu initial state subprocess at the LHC, the black 
lines represent the LO results, while the red ones represent the NLO results. 

Figure UM shows the transverse momentum distributions of the leading jet for the single 
top production via the FCNC couplings together with the SM t-channel single top produc- 
tion, which is the main SM background of our process. We can see that the distributions of 
the single top production via the FCNC couplings decrease more quickly than the SM ones 
with the increasing pt, and the NLO corrections increase the distributions of the FCNC 
single top production, especially in the low pt regions. In Fig. (TT] we show the pseudorapid- 
ity distributions of the leading jet. For the FCNC single top production, the distributions 
decrease with the increasing pseudorapidity at both the Tevatron and LHC, while for the 
SM ones the distributions are almost flat at the Tevatron and increase at the LHC. The 
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NLO corrections increase the distributions by almost the same amount in all the regions. 




p Tj (GeV) P Tii (GeV) 

FIG. 10: Transverse momentum distributions of the leading jet, the black and red lines represent 
the LO and NLO results of the FCNC single top production, respectively, while the green lines 
correspond to the SM ones normalized to arbitrary units. 



Tevatron 



teg 
tug 



1 ' : 


LHC 


teg 
tug 




, ' 1 


i |^ 























InJ 



FIG. 11: Pseudorapidity distributions of the leading jet, the black and red lines represent the LO 
and NLO results of the FCNC single top production, respectively, while the green lines correspond 
to the SM ones normalized to arbitrary units. 

Fig. [12] gives the top quark energy distributions. We can see that the SM ones decrease 
faster than the one of tug FCNC coupling with the increasing top quark energy. Figure [T3l 
shows the invariant mass distributions of the top quark and the leading jet. The shapes for 
the FCNC single top production are different from the SM ones, where there is a peak in the 
middle region. The NLO corrections do not change the shapes of these two distributions. 
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FIG. 12: Energy distributions of the top quark, the black and red lines represent the LO and NLO 
results of the FCNC single top production, respectively, while the green lines correspond to the 
SM ones normalized to arbitrary units. 
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FIG. 13: Invariant mass distributions of the leading jet and the top quark, the black and red lines 
represent the LO and NLO results of the FCNC single top production, respectively, while the green 
lines correspond to the SM ones normalized to arbitrary units. 



Figure [H] shows the jet multiplicities distributions of the FCNC single top production. 
At the LO there is only one jet in the final state, while at the NLO there may be two jets. 
For example, at the LHC, for the tug coupling the LO one jet cross section is 6.7 pb, and the 
NLO corrections reduce the one jet cross section to 4.8 pb and increase the two jets cross 
section to 3.7 pb, which is about 50% of the LO total cross section. 
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FIG. 14: Jet multiplicities distributions of the FCNC single top production, the black and red lines 
represent the LO and NLO results, respectively. 



V. CONCLUSIONS 



In conclusion, we have investigated the NLO QCD effects on the single top productions 
induced by model-independent tqg FCNC couplings at both the Tevatron and LHC. Our re- 
sults show that, for the teg coupling the NLO corrections can enhance the total cross sections 
by about 60% and 30%, and for the tug coupling by about 50% and 20% at the Tevatron 
and LHC, respectively, which means that the NLO corrections can increase the experimental 
sensitivity to the FCNC couplings by about 10%— 30%. Moreover, the NLO corrections re- 
duce the dependence of the total cross sections on the renormalization or factorization scale 
significantly, which lead to increased confidence on the theoretical predictions. Besides, we 
also evaluate the NLO corrections to several important kinematic distributions, i.e., the 
transverse momentum and pseudorapidity of the leading jet, energy of the top quark, jet 
multiplicities, and invariant mass of the leading jet and top quark. We find that for most 
of them the NLO corrections are almost the same and do not change the shape of the 
distributions. 
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